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^ \ Abstract 

^ ] We present the results of a numerical calculation of the B — ► K*^ 

form factors. The results have been obtained by studying the relevant 
correlation functions at (3 = 6.0, on an 18 3 x 64 lattice, using the 0(a)- 
improved fermion action, in the quenched approximation. Prom the study 
of the matrix element {K*\sa flu b\B) we have obtained the form factor Ti(0) 
which controls the exclusive decay rate. The results are compared with the 
recent results from CLEO. We also discuss the compatibility between the 
scaling laws predicted by the Heavy Quark Effective Theory (HQET) and 
pole dominance, by studying the mass- and g 2 -dependence of the form 
factors. From our analysis, it appears that the form factors follow a mass 
behaviour compatible with the predictions of the HQET and that the (in- 
dependence of T2 is weaker than would be predicted by pole dominance. 
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Rome, Italy. 

2 Laboratoire associe au Centre National de la Recherche Scientifique. 



1 Introduction 



An important class of 5-meson decays is given by the weak radiative decays B — > 
Xg'-f, where X s is a strange hadronic state and where the emitted photon is real. 
These decays are extensively studied because they provide, through loop effects, 
interesting information on the Standard Model parameters (e.g. the combination 
of CKM matrix elements |KsV^|). Short-distance physics puts in the foreground 
an effective magnetic interaction, b —>■ s'j, originating from the so-called penguin 
diagrams. In these diagrams, the top quark dominates, whereas the charm and 
up quark contributions are suppressed by powers of the quark masses. 

Radiative decays are also sensitive to physics beyond the Standard Model, 
through possible extra particles (SUSY particles, extra charged Higgs bosons, ...) 
contributing to the loops. Should there be deviations between the expected decay 
rates and the measured values, these could be a manifestation of new physics. 

Among radiative decays, the process B — > K*^ has received an increasing 
attention, because of the experimental measurement of the B — > K*^ branching 
ratio, performed by the CLEO collaboration 0: BR(B — > K*j) = (4.5 ± 1.5 ± 
0.9) x 10 -5 . More recently, the same collaboration has also measured the inclusive 
rate, obtaining BR(B -> X sl ) = (2.32 ± 0.51 ± 0.29 ± 0.32) x 10" 4 @. 

Several methods have been employed to predict the inclusive B — > X s ^ and 
exclusive B — > K*j decay rates: Heavy Quark Effective Theory (HQET)||, QCD 
sum rules [H-|§, and quark models 0-||n]|. For the exclusive decay, the theo- 
retical uncertainty, which was originally of more than two orders of magnitude, 
has been greatly reduced by the more recent studies. Still, there is a large spread 
between the different results. Lattice QCD offers the possibility to investigate 
rare B decays from first principles. The feasibility of the lattice approach was 
first demonstrated in 1991, by the work of Bernard, Hsieh and Soni |12| . Further 
results have been obtained by the UKQCD collaboration [T3j and by the LNAL 
collaboration HI • 



In this paper, we present a lattice calculation of the form factors T\ and T 2 
relevant for exclusive B — ► K*^y decays, and of the dependence of these form 
factors on the heavy-quark mass and on the squared momentum transfer q 2 . 
Preliminary results of our study can be found in ref. jl5fl . The study of the 
dependence of the form factors on the heavy-quark mass provides a good test 
of the validity of the scaling laws predicted by HQET, in a region of masses 
around the charm quark mass. Particular attention has been devoted to the 
understanding of these scaling laws and their relations to the q 2 - dependence. 
This point will be discussed in detail in section |. 

In the same context, there is a very interesting challenge: assuming HQET and 
SU(3) symmetry, the hadronic matrix elements (K*(rj, k)\sa^ l/ q v b\B{j))) are re- 
lated to those relevant for 5-meson semileptonic decays: {p(rj, k)\sj^b\B(p)) and 
(p(r},k)\s r y tlr Y5b\B(p)) |T^|-[T^|. This implies relationships between the semilep- 
tonic and radiative form factors. These relations can be proved, in the infinite 
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mass limit, near the zero recoil point, i.e. q 2 ~ q^ax- A complete test of these 
relationships on the lattice is a very interesting check of QCD dynamics. 

With some extra assumptions one can extend the HQET relations, to the 
region of small q 2 ||18|-||20|. As an example, consider the ratio 



R(B^K*-i) 192tt 3 1 (m 2 B -m 2 K : 



dr(B ^f Pe) \ q z=0 " G F \ V ub\ 2 K-m^ ( m 2- m 2)3 



A' , (1) 



where 



P (R , r(£? — » K*j) (m B + m p ) 2^(0) 

and <^ p (0) is one of the six semileptonic form factors (see for example 
The quantity X can be measured on the lattice. According to [18fl— 1|20||, X should 



be close to 1. Unfortunately, the quality of our numerical results does not allow 
a calculation of X at present. For this reason, we only focus here on B — > K*j 
decay. 

We present results obtained by assuming different scaling laws, in the heavy- 
meson mass, for the form factors at q 2 = 0. These scaling laws will be discussed 
in detail below. The two possibilities correspond to a pole-dominance behaviour 
in q 2 either for T\ or for T 2 and lead to quite different results. 

If we assume that T\ follows the pole-dominance behaviour, we expect a scal- 
ing law of the form Ti(q 2 = 0) = T 2 (g 2 = 0) ~ m -1 / 2 , where m is the mass of the 
heavy quark. In this case, following ref. [EH]], we obtain 

Ti(g 2 = 0) = 0.23 ±0.02 ±0.02, 
R = 0.31 ±0.12, 
BR(B -> K*i) = (7.4 ± 1.4+ 2 ;^) x 10~ 5 . (3) 

where R = T(B — > K* / y)/T(B — > X a j). The formulae used to compute R 
and BR(B —>■ K*j) can be found in section j|. If instead, we assume that 
T 2 follows the pole dominance behaviour, we expect a scaling law of the form 
T\(g 2 = 0) = T 2 (q 2 = 0) ~ m" 3 / 2 . In this case we find 

T 1 (q 2 = 0) = 0.09 ±0.01 ±0.01, 
R = 0.05 ±0.02, 
BR^B-^K^) = (1.1 d= 0.3±g:f) X 10~ 5 . (4) 

The values of Ti(g 2 = 0) in eqs. (|3|) and (|]) are in good agreement with the 
results of other similar studies jl2|-|l4l. In the first case, eqs. (Q), the agreement 
of the lattice predictions with the experimental measurements is rather satisfac- 
tory. In the second case, eqs. (§), either there is a problem with the lattice 
calculations, or the difference may come from long-distance contributions to the 
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exclusive decay rate J22|. Unfortunately, the present lattice technology cannot 
estimate these (and others) long-distance contributions. Our study of the mass- 
and g 2 -dependence of the form factors favours solution (|^), even though, given 
the statistical accuracy of our results and the systematic errors, we cannot draw 
any firm conclusion at this stage. 

The systematic difference between results (|j) and (|4j) originates from the ex- 
trapolation to large meson masses and small q 2 , which is needed to obtain the 
physical form factors. This problem is intrinsic to the lattice approach at values of 
the lattice spacing currently accessible in numerical simulations. In this respect, 
at present, the lattice approach is not very different from quark- model calcula- 
tions. The q 2 - and mass-dependence of the form factors, including those relevant 
to semileptonic decays, remains a crucial challenge for lattice calculations. 



2 Effective Hamiltonian and notation 

In this section, we introduce the essential notation and express the matrix ele- 
ments of the effective Hamiltonian and the exclusive decay rate in terms of the 
relevant form factors. The operator basis of the effective Hamiltonian density, 
for weak radiative 5-meson decays, consists of local four-quark operators O n 
{n = 1, 2, . . . , 6) and magnetic-type operators O n {n = 7, 8) p6|-pT|: 

H eS = -V tb V*G F /V2j2Cn(fi)O n (fi), (5) 

n=l 

where C n are the Wilson coefficients. The operator that controls the b —>■ s'j 
transition is: 

e 

•16V 2 

From eqs. (|5|) and (j^), one finds the matrix element for the transition B — > K*^ 
M = e ^0C 7 (m b )V tb V t ^(K*\J,\B) , (7) 

where 

and & and q^ = — k^ are the photon polarization and momentum transfer, 
respectively. We parametrize the hadronic matrix element in eq. (|7|) as 

(K*(r h .,k)\J»\B(p)) = C?T 1 (q 2 )+ l CZT 2 (q 2 )+ l C»T 3 (q 2 ), (9) 



7 =(—)m b (7 S( r^b) R F llu . (6) 



where 



CT = 2e^rj* r (k) aPp k a , 

C% = rjrikXM'-M^-^ik^ip + ky, (10) 
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Ti, T 2 and T3 are real dimensionless Lorentz-invariant form factors, and 77 is the 
K* polarization vector. The vector current Jy contributes only to Ti, the axial 
current J'X only to T 2 and T 3 : 



(K*( Vr ,k)\Jy\B( P )) 

(K*( Vr ,k)\r A \B(p)) 



iC^T 2 {q 2 ) + iC^T,{q 2 ) 



T 3 does not contribute to the physical rate, because its coefficient vanishes for 
a transversely polarized photon. By using the relation a M ^7 5 = \e^ ual3 a a p, one 
finds 



T 1 (0)=T 2 (0) 

at q 2 = 0. From the matrix element (|7|), using ( 



(12) 

one obtains the decay width 



F(B -> iT 7 ) = j^rnlGlMl f 1 



M 



K* 



Ml 



VaVffCrimtflT^l 2 . (13) 



The physics of this decay is then described by one form factor only, T\ at q 2 = 0. 

It is not convenient to compare the theoretical width ( |T3"D with its experi- 
mental value, because many theoretical uncertainties enter in this quantity: the 
renormalization scale, the matching of the lattice to the continuum operators, the 
value of Aqcd, the possible presence of new physics beyond the standard model, 
uncomputed higher-order corrections, etc. |21|| . Now that the inclusive measure 
is available ||, it is much more informative to compare instead the exclusive- 
to-inclusive ratio of rates R because in this ratio most of the above-mentioned 
uncertainties cancel out. In terms of Ti, the ratio R can be written as 



R 



T(B - X sl ) 



th 



Mb 



1 - 



Ml 



x 



1 + (Ai - 9A 2 )/(2m£ 



•xlTx 



(14) 

where we have divided the exclusive rate flT3|) by the theoretical inclusive rate, 
computed in the HQET parton model. In this formalism, the parameters Ai and 
A 2 describe the leading non-perturbative corrections (at order 0(1 /m 2 )) to the 
parton-model predictions for the inclusive rateQ f32fl . 

The branching ratio BR(B — > K*j) is conveniently expressed through a chain 



of ratios [21 



BR(B -> K*j) = R x 



T(B -> X sl ) 



tii 



x BR exp {B -> Xlvi). 



(15) 



1 They are related to the kinetic energy of the &-quark inside the B-meson and to the B-B* 
mass splitting. 
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3 Correlation functions and extraction of the 
form factors 



In lattice simulations, from the study of two- and three-point correlation func- 
tions, one extracts the current matrix elements for a given momentum transfer 
and for a given polarization rj of the K* . This is by now a well established tech- 



nique. The reader can find more details in refs. |33| and ||34|| . We follow more 



closely ref. [33 



In general, the form factors are obtained at values of q 2 that are constrained 
by the quark masses and the values of particle momenta accessible on a given 
lattice volume. For this reason, it is necessary to extrapolate T\ and T 2 to q 2 = 0, 
in order to get the physical result. In the following, the different steps of the 
procedure to extract the form factors are briefly described^. 

The matrix elements have been computed for an initial meson at rest and a 
final vector meson with momentum pk*- We have taken p^* = 2ii/(La) (0, 0, 0), 
2vr/(La)(l, 0, 0), 2vr/(La)(l, 1, 0), 2vr/(La)(l, 1, 1), and 2vr/(La)(2, 0, 0); where L 
is the spatial extension of the lattice, in our case L = 18, and a is the lat- 
tice spacing. Correlation functions, which are equivalent under the hypercubic 
symmetry, have been averaged. 

The initial (final) meson was created (annihilated) by using the pseudoscalar 
(Jb = bj^q) and local vector (J^* = q^ a s) densities, inserted at times t^/a = 28 
and tx* = 0, respectively. We have varied the time position of J M in the interval 
tj/a = 10-14. 

In order to obtain the hadronic matrix element, the following procedure has 
been used: 

1. We have computed masses and source matrix elements by fitting the two- 
point correlation functions at large time distances to the expressions given 
below. For the pseudoscalar meson, we have used: 

C B (t B ) = (0|J B (t B )4(x = 0,t = 0)|0) = ^e- M ^, (16) 

where Jb{^b) = J d 3 x Jb{x, t#) and Z^ 2 = (0\Jb\B(pb = 0)). For the 
vector current, one has 



where 



C#(p K ;t K .) = (0\J%.(pK;tK')f&@=0,t= 0)|0) 

( a/9 _ PK* A' \ Z K* —Ere* t K * /i - 

m 2 k , ' 2Ek* ' 1 

J%*{p K *,t K *) = Jd 3 xe- ip ~«* s J%*(x,t K *). (It 



2 We adopt the convention that B and K* denote the pseudoscalar heavy meson and the 
vector light meson, whenever no confusion arises. 
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In eq. (PI), Zk* is denned through 



r^y/Z K . = (K*( Vr ,p K *)\J K *\0) . (19) 

2. We have then extracted the matrix elements from the ratios 

P Cg a (t Bl tj,p K *) 



with 



Ct(t B ,tj,p K *) = (0\J K *(p K *,0)J»(-p K *,tj)J B (t B )\0) 

V n a ^K* r -EK*t., VZb M B (t B -t. ; ) 

2 ? Vr 2E K * e 2M B 



x (K*(7ir,p K .)\J"\B) , (21) 

where the magnetic operator J 11 , renormalized in the continuum, is given in terms 
of the corresponding lattice bare operator by J M = Z a Jl^ ttl with Z a = 0.98 taken 
from perturbation theory p5 |. To obtain the matrix elements, we have used 



two different procedures, denoted by "analytic" and "ratio" methods, which are 



discussed in detail in ref. |]33| . In the "ratio" method, for each fixed-time distance, 
the three-point correlation function is divided by the two-point functions of the 
B and K* mesons (with corresponding momentum) as in eq. (|20|), in order to 
cancel the exponential time-dependence. The "analytic" method differs from the 
previous one in that, instead of dividing by the numerical two-point correlation 
functions, we divide by the corresponding analytic expressions using the source 
matrix elements (Z B and Zk*) and the meson energies obtained from the fit of 
the two-point functions at zero momentum. The energy of the K* si computed 
from the meson mass, as explained in eqs. (^) and (|25|) below. When both the 
initial and final mesons are at rest, and at large time distances, the two methods 
are practically equivalent and lead to almost identical results. However, when the 
meson momenta are different from zero the two methods are expected to agree 
only up to 0(a) effects. 



In refs. |T3, Rl, RBI, it was found that discretization errors appear to be re- 



duced, if one uses, for the two-point correlation functions, the "lattice" dispersion 
relation of a free boson 

= (22) 

where 

E = -arcsinh I sinh 2 ( ^ ) + V sin 2 ( ^ ) ] (23) 




and p is the momentum of the meson. The same is true in our case. We have 
verified this point by studying the ratio 7Z(t) = C(t,p)/C(t,p), where C(t,p) is 
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the two-point correlation function of a meson with momentum p, as computed 
in our simulation, and C(t,p) is given by the expression in eq. (|2~2"D , with Z and 
m taken from the fit of the zero- momentum correlation to C(t,p = 0). At large 
time distances, we find that 7l(t) is closer to 1 (typically 1.05 ± 0.01 instead of 
1.10 ± 0.02), if we use eq. ( p2"D instead of the standard expression 

C(t,p) = ^e- E \ (24) 

with 

E = ym 2 + \p\ 2 . (25) 

For this reason, throughout our analysis, we have fitted the two-point functions 
to C(t,p), eq. (p2|), in addition to the "standard" form C(t,p), eq. fl2"4"|). 

Using C(t,p), the "ratio" and "analytic" methods yield only slightly different 
results (see section ||). 



4 ^-dependence of the form factors and scaling 
laws 

In order to obtain the form factors at the physical point, we need to extrapolate in 
both mass and momentum. The final results depend critically on the assumptions 
on the q 2 - and heavy mass-dependence, which deserve a detailed discussion. 
At fixed Pk*-, with \pk*\ *C Mb in the 5-meson rest frame, the following 



scaling laws can be derived fL6f : 



these are valid up to logarithmic corrections. As mentioned in the introduction, 
"scaling" laws for the form factors at q 2 = can only be found by using ex- 
tra assumptions on the (^-dependence of the form factors. Such a procedure is 
acceptable, provided the "scaling" laws derived in this way respect the exact con- 
dition Ti(0) = T2(0). This is a non-trivial constraint: the g 2 -behaviour of 7\ and 
T2 has to compensate the different mass dependence of the two form factors near 
the zero recoil point, see eq. (^B|). Thus, for example, the popular assumption of 
pole dominance for both Ti and T 2 would give 71(0) ~ Mg 1 ^ 2 and T 2 (0) ~ M B 3 ^ 2 , 
which is inconsistent with eq. (p~2|) . 

In all lattice simulations, which try to compute i?-meson form factors by 
extrapolating from lower heavy-quark masses, we are forced to make assumptions 
on the corresponding g 2 -dependence. This is a consequence of the fact that the 
extrapolation in the mass, at fixed light-meson momentum, pushes the values of q 2 
towards q max . This problem can only be avoided by going to much smaller values 
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Figure 1: The inverse form factor 1 / (2T2(g 2 )) is shown as a function of q 2 /M 2 for 
Kh = 0.1200. The form factor has already been extrapolated in the mass of the 
light quark to the strange-quark mass. The dash-dotted line represents the pole- 
dominance behaviour, with the mass of the axial meson M t taken from the fit of 
the two-point correlation function; the solid line corresponds to a pole- dominance 
behaviour, with M t left as a free parameter; the dashed line is a fit with T2(q 2 ) 
constant in q 2 . 



of the lattice spacing, thus allowing a large increase in the range of accessible 
masses and momenta. The assumptions on the g 2 -dependence of the form factors 
can be tested, although in a small domain of momenta, directly on the numerical 
results. 

In fig. |l|, we show l/(2T 2 (g 2 )) as a function of the dimensionless variable 
q 2 /M 2 , where M t is the appropriate mass for the axial meson exchanged in the 
t-channel, in the pole-dominance approximation. If pole dominance is valid, the 
form factor should have a linear behaviour in q 2 , with a slope proportional to 
1/M 2 . In the case of T 2 , we found that the mass extracted from the slope is 
much larger (by a factor of order 2 in most of the cases) than the mass that 
we have obtained from the axial two-point correlation function. In other words, 
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T 2 (g 2 ) is flatter than predicted by pole dominance. In the case of T±, the point at 
momentum (0, 0, 0) cannot be computed. The absence of this point and the large 
statistical and systematic errors at large momenta make it difficult to test the 
validity of the pole- dominance hypothesis directly on our results. A consistency 
check will be provided is section [5]. Although our data are not accurate enough 
to draw a definite conclusion, they suggest that assuming T 2 almost constant in 
q 2 and T\ following pole dominance gives a good description of our results. This 
assumption is consistent with the "scaling" laws described above. We call this 
option m -1 / 2 -scaling. The m~ 1//2 -scaling is similar to what has been encountered 
in the QCD sum-rules calculation [ 37|, |38| of the semileptonic B — > p form factors 
V and Ai, which are related to T\ and T 2 in the large-mass limit. In fig. [I], we 
also give the curves corresponding to a pole dependence for T 2 (m~ 3 / 2 -scaling). 
This possibility is also consistent with the "scaling" laws, if a dipole dominance 
behaviour is assumed for 7\, but it is disfavoured by our numerical results, as 
discussed above and shown in the figure. For these data, the uncorrelated y 2 , 
corresponding to the fit of T 2 to a constant, is x 2 = 0.65 to be compared to 
X 2 = 1.15, in the case of the fit of T 2 to a pole dominance behaviour. Similar 
results have been obtained at the other values of the heavy quark mass. It 
is clear, that an infinite number of possible g 2 -dependence for 7\ i2 , which are 
compatible with the scaling laws, can be found. However it would be impossible 
to distinguish among them, given the statistical errors, the systematic uncertainty 
in the extraction of the form factors, the effects of 0(a), and the limited range in 
q 2 and Mb- Thus we take the two possibilities, mT 1 ' 2 -scaling and m~ 3//2 -scaling, 
as representative of a full class of "scaling" laws. In section |5|, we will see that the 
two options lead to quite different results for the physical value of Ti(0), which 
enters in the calculation of the B — > K*^ rate. 

There remains to discuss the dependence on the light-quark masses (spectator 
and active). At fixed heavy-quark mass and light-meson momentum pk*, the 
generic form factor F (F = T 1; T 2 ) has been extrapolated linearly in the mass 
of the active light quark 

F = a + pm q , (27) 

to values corresponding to the physical strange meson K*, assuming the form 
factors independent of the mass of the spectator quark. This is in agreement 
with the results of ref . [F3 , where it was shown that the dependence of the form 
factors on the mass of the spectator quark is very small. We believe that the 
extrapolation in the light-quark mass is unlikely to be a source of an important 
uncertainty within the present statistical accuracy. 



5 Lattice setup and numerical results 

The numerical simulation was performed on the 6.4 gigaflops APE machine, at 
(3 = 6.0, on an 18 3 x 64 lattice, using the SW-Clover action in the quenched 
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Form factor 


Kh 


r-sinh 


a-sinh 


r-stand 


a- stand 


Ti 


0.1330 


0.35(3) 


0.34(3) 


0.33(3) 


0.34(3) 


T 2 


0.1330 


0.37(2) 


0.37(2) 


0.35(2) 


0.35(2) 


f 2 


0.1330 


0.37(2) 


0.36(2) 


0.35(2) 


0.35(2) 


T x 


0.1250 


0.34(4) 


0.33(3) 


0.31(3) 


0.32(3) 


T 2 


0.1250 


0.36(2) 


0.36(2) 


0.33(2) 


0.33(2) 


T 2 


0.1250 


0.34(3) 


0.33(3) 


0.31(3) 


0.32(3) 


T x 


0.1200 


0.33(4) 


0.32(3) 


0.30(3) 


0.31(3) 


T 2 


0.1200 


0.35(2) 


0.35(2) 


0.32(2) 


0.32(2) 


T 2 


0.1200 


0.33(5) 


0.32(4) 


0.29(4) 


0.31(4) 


Ti 


0.1150 


0.33(4) 


0.32(4) 


0.29(3) 


0.30(3) 


T 2 


0.1150 


0.35(2) 


0.35(2) 


0.30(2) 


0.30(2) 


T 2 


0.1150 


0.32(6) 


0.31(6) 


0.28(5) 


0.30(6) 


Ti 


B 


0.24(3) 


0.23(3) 


0.20(3) 


0.21(3) 


T 2 


B 


0.25(2) 


0.25(2) 


0.22(2) 


0.22(2) 


T 2 


B 


0.22(5) 


0.21(5) 


0.19(5) 


0.21(5) 



Table 1: T 1; T 2 and T 2 at q 2 = for different values of the heavy-quark masses, 
as extracted using the pole and constant behaviour respectively. We give in the 
second column the hopping parameter of the heavy quark (B denotes the extrapo- 
lation to the B-meson using the mr 1 ! 2 scaling law); the light-quark mass has been 
extrapolated to the strange quark mass; the values of the form factors obtained 
with the ratio method from the sinh or standard fits (r-sinh and r-stand) or with 
the analytic method from the sinh or standard fits (a-sinh and a- stand) are shown 
from the third to the sixth columns. 
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0.29(3) 


T 2 


0.1200 


0.29(2) 


0.29(2) 


0.26(2) 


0.27(2) 


Ti 


0.1150 


0.28(4) 


0.28(3) 


0.25(3) 


0.26(3) 


T 2 


0.1150 


0.27(2) 


0.27(2) 


0.24(2) 


0.24(2) 


Ti 


B 


0.10(1) 


0.09(1) 


0.08(1) 


0.09(1) 


T 2 


B 


0.09(1) 


0.09(1) 


0.08(1) 


0.08(1) 



Table 2: T\ and T 2 at q 2 = for different values of the heavy-quark masses, as 
extracted using the dipole and pole behaviour respectively. We give in the second 
column the hopping parameter of the heavy quark (B denotes the extrapolation 
to the B-meson using the m~ 3 / 2 scaling law); the light quark mass has been ex- 
trapolated to the strange-quark mass; the values of the form factors obtained with 
the ratio method from the sinh or standard fits (r-sinh and r-stand) or with the 
analytic method from the sinh or standard fits (a-sinh and a-stand) are shown 
from the third to the sixth column. 



T(0) = T 1 (g 2 = 0) = T 2 (g 2 = 0) 


Scaling law 


r-sinh 


a-sinh 


r-stand 


a-stand 


m -l/2 
m -3/2 


0.25(2) 
0.09(1) 


0.25(2) 
0.09(1) 


0.21(2) 
0.08(1) 


0.22(2) 
0.08(1) 



Table 3: T(0) = Ti(q 2 = 0) = T 2 (q 2 = 0) extrapolated to the B-meson using the 
vrT 1 ! 2 and the m~ 3 ^ 2 scaling laws; the values of the form factors obtained with 
the ratio method from the sinh or standard fits (r-sinh and r-stand) or with the 
analytic method from the sinh or standard fits (a-sinh and a-stand) are shown 
from the second to the fifth columns. 
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approximation. The results were obtained from a sample of 170 gauge config- 
urations. The statistical errors have been estimated by a jacknife method, by 
decimating 10 configurations from the total set. For each configuration we have 
computed the quark propagators for 7 different values of the Wilson hopping 
parameter Kw, corresponding to "heavy" quarks, Kh = 0.1150, 0.1200, 0.1250, 
0.1330, and "light" quarks, K L = 0.1425, 0.1432 and 0.1440. In order to extract 
masses and source matrix elements, we have fitted the heavy and light meson 
two-point functions to eqs. (|16|) and (|TT|) in the time interval 14 < t/a < 28 
and 10 < t/a < 28 respectively. The relevant matrix elements of O7 have been 
extracted, for different polarizations and momenta of the K*, from the three- 
point function, see eq. (|2l|) , with the operator inserted at 10 < tj/a < 14. 
We found that the critical value of K^, corresponding to the chiral limit, is 
K cr = 0.14545(1); the inverse lattice spacing, obtained from m p , is a -1 = 1.96(7) 
GeV; the value of the Wilson parameter for the strange quark, determined from 
the mass of the kaon, is K s = 0.1435(1). In the following, the labels "ratio" and 
"analytic" refer respectively to the ratio and analytic method as explained in 
section |3|. The numbers quoted with "sinh" and "standard" are obtained respec- 
tively with the fits to eqs. ( p2|) and (|24]) for the two-point functions (see section 

We now explain the analysis of the scaling behaviour of the form factors that, 
on the basis of the above discussion, has been done in combination with the 
study of their g 2 -dependence. Unless otherwise stated, the form factors are those 
obtained after the extrapolation in the active light quark mass to the value of the 
hopping parameter corresponding to the strange quark. 

1. According to HQET, when the mass of the meson is sufficiently large, T 2 
at zero recoil should follow the behaviour given in eq. (|26|). A fit of our 
data for T 2 (q^ nax ) to MJj(a + b/M B ), with a, a and b as free parameters, 
gives a = -0.41(10), and b/a = -350(70) MeV (a = -0.68(10), and 
b/a = —480(50) MeV) with the sinh-fit (standard- fit). The exponent a 
is thus compatible with the value of —1/2 predicted by HQET, see eq. 
(|26|). To reduce the number of parameters, we have then extrapolated 
T2((lm ax ) in 1/Mb, with a = —1/2, using eq. (|2"6]). In this case, we obtain 
72 = 20.2 (1.6) MeV 1 / 2 and S 2 = -430 (50) MeV (72 = 16.8 (1.4) MeV 1 / 2 
and 82 = —320(50) MeV). These parameters correspond to the following 
values of T 2 ( y q^ nax ) for the B- meson 

sinh standard 
T 2 (q 2 max ,B)= 0.25(2) 0.22(2). (28) 

This result is also reported in table [IJ as T 2 with the label B. 



12 



2. Let us assume the "pure-m -1 / 2 " behaviour for the form factors. By "pure- 
m -i/2» behaviour, we mean that T 2 (q 2 ) is independent of q 2 (this is com- 
patible with our results, cf. fig. [1]) and that Ti(q 2 ) follows a pole-dominance 
behaviour, with the mass of the vector meson as measured from the two- 
point correlation functions. In table |1], we compare, for all the values of 
the heavy-quark mass, T 2 = T 2 (q 2 = 0) = T 2 (q^ nax ) with Ti = Ti(g 2 = 0) = 
T\{q 2 ) x (1 — q 2 /M 2 ). The values of Ti and T 2 are compatible within the 
statistical errors and the differences coming from different extrapolations 
(analytic, ratio, sinh, standard). 

3. In order to check the stability of the results, we also performed a linear 
extrapolation in q 2 of T 2 (q 2 ) to q 2 = 0, at fixed heavy-quark mass, followed 
by a fit of T 2 (q 2 = 0) to eq. fl2"6|). The values obtained in this way are 
reported as T 2 in table p]. Since we have not done any assumption on 
the Me-dependence of the slope of the fit, this procedure is not a priori 
incompatible with the scaling laws discussed in the previous section. It 
should be noticed that, since we have an extra parameter in the fits, i.e. 
the linear slope in q 2 , the extrapolated value T 2 has a larger error. 

4. Even though our results prefer a flat g 2 -dependence for T 2 , we have also 
fitted Ti and T 2 by assuming a pure "m -3 / 2 " behaviour for the form fac- 
tors. This means that we first fit the g 2 -dependence of the form factors to 
Ti(g 2 ) = Ti/(1 - q 2 /M' t 2 ) 2 and T 2 (q 2 ) = T 2 /(l - q 2 /M 2 ), where for each 
value of the heavy-quark mass, the value of M[ and M t are those computed 
from the corresponding two-point functions. The form factors extracted in 
this way are reported in tab. |2|. We have then extrapolated in 1/Mg using 
the expressions Ti M 3 B /2 = %i x (l+0i/Af B ) and T 2 M 3 J 2 = X 2 x {l+6 2 /M B ). 

5. We notice that, for the values of the heavy quark masses at which we have 
computed the form factors, the value of Ti(g 2 = 0) (T 2 (q 2 = 0)) extracted 
by assuming the pole (constant) g 2 -behaviour differs at most by 20-25% 
from the value obtained by assuming the dipole- (pole-) g 2 -depencence, as 
can be read in tables [l] and |2|. The same would be true if we extrapolated 
the form factors using the two different scaling laws to the charm-quark 
mass. The values extrapolated to the 5-meson instead, differ by about a 
factor of 2. 

6. One should also take into account the distortion in the mass dependence of 
the form factors due to lattice artefacts at large quark masses. This effect 
has been measured non-perturbatively both by the APE collaboration at 
13 = 6.0 gg and the UKQCD collaboration at = 6.2 jjlj. Only with more 



accurate data, and exploring a much larger range of masses and momenta 
at larger values of j3, will it be possible to decide the fundamental issue of 
the scaling behaviour. 
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7. We can exploit the information that the two form factors must be equal 
at q 2 = 0, by making a combined fit of T\ and T 2 , with the constraint 
T(0) = Ti(q 2 = 0) = T 2 (g 2 = 0). This can be done both in the m _1//2 and 
m~ 3 / 2 cases. The results are given in table [| 

8. As a consistency check of our results, we have also inverted the order of the 
fits, by fitting first T 2 (q^ aax ) in 1/M B and then assuming the q 2 behaviour 
of the form factor. The results are, within the errors, well compatible with 
those reported in tables [I] and 0. 

From the results reported in tables we quote 

T!(g 2 = 0) = 0.23(2)(2) scaling m~ 1/2 , 

Tx(g 2 = 0) = 0.09(1)(1) scaling m" 3/2 , (29) 

from which we have derived the results of (|3[) and (|j) in the introduction. 



6 Conclusion 

We have computed the form factor relevant for B — > K*^ decays. In order to 
extract the physical form factor from the lattice results, we have extrapolated 
in momentum transfer and in the mass of the heavy quark. The results, corre- 
sponding to different choices of the scaling law, can differ by about a factor of 
2. Within large uncertainties, our study suggests that the scaling law "m~ 1//2 ", 
by which T 2 (g 2 ) has a very small dependence on the momentum transfer, is pre- 



ferred, in agreement with the results of ref. [[L4J]. We cannot exclude, however, 
the "m~ 3//2 " scaling behaviour, or any intermediate solution. In order to improve 
the accuracy of the predictions, it is necessary to be able to work with heavier 
quark masses and to increase the range of q 2 . This can only be achieved by going 
to larger lattices. 
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